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Abstract: The theory accounting for nonlinear particle relaxation (viscous and thermal) has been 
applied to the prediction of dispersion of sound in dilute suspensions. The results suggest that 
significant deviations exist for sound dispersion between the linear and nonlinear theories at large 
values of cot 4 , where <51 is the circular frequency, and r^ is the Stokesian particle relaxation 

time. It is revealed that the nonlinear effect on the dispersion coefficient due to viscous 
contribution is larger relative to that of thermal conduction. 

PACS numbers: 43.20.Hq, 43.50.Nm, 43.50.Gf 

1. Introduction 

Sound attenuation in fluids, representing the dissipation of acoustic energy from a sound wave, 
occurs through a number of physical processes involving molecular viscosity, thermal 
conductivity, and other dissipative or relaxation processes. 1 ' 7 When a fluid contains 
inhomogeneities such as suspended particles (solid particles, drops and bubbles) additional 
viscous and heat conduction losses occur in the immediate neighbourhood of the suspended 
particles. 1 * 3,8 Particle relaxation in a dilute suspension is the process by which particles adjust to 
fluctuations (velocity and temperature) of the surrounding fluid, leading to attenuation and 
dispersion of a sound wave. A comprehensive review of the physics and scientific history of 



acoustic interactions with particulate mixtures is provided by Challis et al. 9 The acoustic intensity 
1 of a plane wave propagating through an absorbing medium is expressed by 

/ = v ¥ (1) 

l 

where x is the distance traversed, / 0 the intensity at * = 0 , and a, the intensity attenuation 
coefficient for the medium. The quantity a , depends on viscosity, thermal conductivity, and other 
factors such as molecular relaxation. 

Sound propagation in aerosols and fog has been studied experimentally and theoretically 
by several investigators., since the pioneering work of Sewell 10 , with the aid of a scattering 
formulation on the assumption of immovable particles. Epstein" extended this theory for 
particles in motion, and Epstein and Carhart 12 additionally considered heat conduction effects. 
Allegra and Hawley 13 provided further extensions by including liquid-liquid as well as liquid- 
solid systems. 

Temkin and Dobbins, 14 in their classical work involving a coupled-phase formulation, 
theoretically considered particle attenuation and dispersion of sound in a manner which illustrates 
explicitly the relaxation character of the problem. Basset history and added mass terms were 
included in an elegant coupled phase formulation by Harker and Temple. 15 Coupled phase effects 
were also treated by Evans and Attenborough 16,17 in an extension to the work of Harker and 
Temple 15 to incorporate thermal conduction. 

The absorption of sound in suspensions of irregular (nonspherical) particles was 
considered by Urick 18 . Experimental and theoretical studies (extension of Urick's model) on 
sound absorption involving irregular particles were also considered by Richards et al. 19,20 

The particulate relaxation models for the sound attenuation are all based on Stokes drag 
(linear drag law) and pure conduction limit (linear heat transfer). Recently the author 21 
investigated sound attenuation in dilute suspensions, and extended the theory of Temkin and 
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Dobbins 14 by considering nonlinear drag and heat transfer laws applicable to relatively large- 
sized droplets. The absorption coefficient per unit frequency predicted by the nonlinear theory is 
compared with that indicated by the theory of Temkin and Dobbins 14 in Fig. 1. In Fig. 1, 



a =c 0 a/{C m a)) = f(a>T d ,Pr,c pp /c pg ,r) 

(2a) 

where 

* d =d 2 p P P 7 ( 18 ^) 

(2b) 

and 

C m = n 0 m p ! Pg 

(2c) 


In the above, a is the attenuation per unit frequency per unit mass fraction, a the amplitude 
attenuation coefficient, cothe speed of sound in the gas phase, rothe circular frequency, C m the 
mass concentration, r d is the dynamic relaxation time of the particle (relating to particle-fluid 
velocity lag), rtQ the mean number of particles per unit volume of mixture, and m p the mass of 
one particle, p p the mean particle density, p g the mean density of gas, p g the mean dynamic 
viscosity of gas, and d p the particle diameter. Also the quantity Pr refers to Prandtl number of 
the gas, c pg ihe specific heat of gas, c pp the specific heat of the particle, and /the isentropic 
exponent (specific heat ratio). Note that a = aj 2. The results shown in Fig. 1 correspond to 
Cpp / c pg = 4.17, Pr = 0.71, y = 1 .4 (representative of water droplets in air 14 ). With the aid of this 

nonlinear model, the existence of the spectral peak in the linear absorption coefficient a has been 
demonstrated 21 . 

On the basis of this extension, good agreement was achieved with the recent data of 
Norum 22 for sound attenuation in perfectly expanded supersonic jets containing suspended water 
droplets, which reveal that the linear absorption coefficient displays a spectral peak (Fig. 2). The 
data correspond to hot supersonic jet of air from a convergent-divergent (CD) nozzle operation at 
a jet total temperature T t = 867 K, and a jet exit Mach number Mj= 1.45. The jet Mach number is 
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defined as M j = uj Icj , where the subscript j refers to the nozzle exit conditions. The mass flow 

rate (maximum considered) of water to that of the jet is about 0.85. The angle 0is measured from 
the jet inlet axis. In the data, water is injected at 45 deg. 

Similar spectral peaks in the attenuation coefficient have been observed in the 
measurements by Krothapalli et al. 23 for microjet injection of water into high speed exhaust jets. 
Spectral peaks in the linear absorption coefficient in dilute suspensions were also noticed in the 
experimental data of Richards et al. 19,20 covering a frequency range of 50-150 kHz. 

In the preceding work by the author 21 results were presented only for sound attenuation 
with nonlinear particle relaxation, and the effect of nonlinearity on the dispersion effects was not 
considered. In practical applications such as exhaust jets with injected water droplets, in addition 
to sound attenuation the dispersion effects (phase velocity changes) are also important from the 
standpoint of sound radiation. The present work applies the nonlinear theory for the dispersion of 
sound in dilute suspensions, and compares the linear and nonlinear theories for sound dispersion. 


2. Sound dispersion 


The present analysis for sound dispersion is similar to that proposed by the author 21 for sound 
attenuation, and will be briefly presented as below. Without any loss of generality the dispersion 
of sound for large particle Reynolds numbers with nonlinear particle relaxation may be expressed 
with the aid of Temkin and Dobbins 14 results as follows: 





1 


/ C m = 
m , 2 

\ + arZ 


+(r- 1 ) 


d 1 


Z£P_ 

\ C PSJ 


1 + a) 2 T?\ 


( 3 ) 


In the above the quantity P is the dimensionless dispersion coefficient, c p the specific heat, c 
the actual speed of sound in two-phase medium (phase velocity), and /the isentropic exponent 
(specific-heat ratio). The subscripts g and p respectively denote the gas and the particle. 
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The relaxation times and r,j correspond to those under nonlinear drag conditions 
(generally representative of large-sized particles). Physically the dynamic relaxation time z dx is a 
measure of the time scale in which the particles follow (respond to) the fluctuations in the fluid 
motion. 3 Likewise, the thermal relaxation time r„ is a measure of the thermal response time of 
the particles to follow the fluctuations in the temperature of the fluid. They are related to the 
relaxation times t d and r, by the relations 21 

*d\ = *d¥ l ( R e P \ T t\ = 2 ( Re P ) ( 4a ) 

where p )=C Dl /C D , ^ 2 l Re p> Rr ) = ^ M l (4b) 


with Cdi standing for the nonlinear drag coefficient, and Nu\ for the nonlinear heat transfer 
(Nusselt number). In the above, the quantity xj is given by Eq. (2b), and 


= 


MpCpp ^CppdpPp 


2nd pkg 
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Also the particle Reynolds number Re^ is defined by 
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(4d) 


where u g and u p denote the velocity of the particle and the gas respectively. Also the quantity 
Pr = c pg /i g / k g stands for the Prandtl number of gas, where stands for the thermal 
conductivity of the gas. 

The drag coefficient and the Nusselt number in Eq. (4b) are defined by 


C D =2F p 


{Pg^ d l u 


Nu = h g d p I kg 


( 5 ) 
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where F p is the particle drag force, and h g Xhe gas-droplet convective heat transfer coefficient. 
For rigid particles, the linear droplet drag and heat transfer are respectively obtained from 

F p = 6np g (u p -u g ) (6a) 

Qp=2ndpk g {T p -T g ) (6b) 

which correspond to the zero droplet Reynolds number limit (Re p -» 0 ). According to Temkin 
and Dobbins, 14 Stokes linear drag law can be justified forQ< cut d »1, provided that 

p g / p p «1 and {codp /Sv g J /2 «1. 

The expressions for y / j and y / 2 have been taken as 21 

\ 

, f„ +0.4 (7) 

+ V Re p j 

^ 2 (Re p ,Pr)= Nu x f Nu = 1 + 0.3 Re J 5 Pr 033 (8) 

Eq. (8) is obtained from Ref. 24, and Eq. (23) is based on the Ref. 25. 

The determination of particle Reynolds number required in the evaluation of the functions 
y/\ and ^ 2 > n Eqs. (7) and (8) respectively is exceedingly complex. There exists relatively little 
information on the dependence of particle Reynolds number on the particle characteristics in two- 
phase flows. In this connection the author 21 postulated that the particle Reynolds number depends 
only on the particle relaxation time, and is independent of the particle to fluid density ratio, and 

the following power law relation is proposed on the basis of work of Ref. 26: 

Re p = f{o)T d ) = c((OT d f ( 13 ) 
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The adjustable constant c is determined from a correlation of the theory with the test data. A 
value of c = 10 was found to be satisfactory based on the data of Norum 21 for water droplets in a 
supersonic air jet (Fig. 2). 

3. Results and comparison 

Figure 3 illustrates a comparison of the predictions for the dispersion coefficient between the 
linear and nonlinear theories for particle relaxation. The results are shown for 
c pp / c pg = 4.17,Pr = 0.71, x = 1 .4 (representative of water droplets in air 14 ). The contributions of 

viscosity and the thermal conductivity along with their combined effect on the dispersion 
coefficient are shown for both linear and nonlinear relaxation. It is seen that the nonlinear effects 
become important for cot j > 0.2 (where heat conduction effects become important). With regard 

to the viscous contribution, nonlinearities become manifest for somewhat higher values of eotj 'm 

excess of about 0.4. The results also suggest that the nonlinear effect is more significant in the 
viscous contribution relative to that of thermal conduction. With nonlinear particle relaxation, the 
total dispersion coefficient approaches zero for smaller values of cotj than those in the case of 
linear relaxation. 

The predictions for the dependence of the dispersion coefficient (yielding the phase 
velocity c) on the frequency for various values of the dynamic relaxation time tj is presented in 

Fig. 4. For example, the dynamic relaxation time tj for silica particles in air is about 10' 3 sec for 
a 5 fan particle, and about 10' 1 sec for a 50 pan particle. 8 On the other hand, the thermal 
relaxation time for particles in air is about 5x1 0 4 sec for a 5 jan particle, and about 5xl0' 2 sec 
for a 50 tan particle. For comparison purposes, the molecular and thermal relaxation time scales 
are of the order of 10' 10 sec for gases. The results suggest that as the particle relaxation time 
increases, the dispersion curve is shifted to lower frequencies, as is to be expected. 
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It is believed the nonlinear particle relaxation effects on sound attenuation and dispersion 
will be of interest in the prediction of jet noise reduction by water injection . 27 * 29 

4. Conclusion 

The theory of nonlinear particle relaxation proposed previously for sound attenuation in dilute 
suspensions has been extended to investigate the effect of nonlinearity on sound dispersion. The 
results reveal that significant nonlinear effects are noticed at relatively large particle relaxation 
times. It is also observed that the nonlinear effect on dispersion due to viscous contribution is 
larger relative to that of thermal conduction. 
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Figure Captions 

Fig. 1. (color online) Predicted absorption coefficient with nonlinear particle relaxation 
processes. 

Fig. 2. (color online) Comparison of the predictions for the linear absorption coefficient 
with test data of Norum 22 . 

Fig. 3. (color on line) Comparison of the predictions for the dispersion coefficient 
between the linear and nonlinear theories for particle relaxation. 

Fig. 4. (color online) Variation of dispersion coefficient with frequency for various values 
of particle relaxation time. 
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Fig. 1. (Color online) Predicted absorption coefficient with nonlinear particle relaxation 
processes. 
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Fig. 2. (Color online) Comparison of the predictions for the linear absorption coefficient with 
test data of Norum 22 . 
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Fig. 3. (Color online) Comparison of the predictions for the dispersion coefficient between 
the linear and nonlinear theories for particle relaxation. 



Fig. 4. Variation of dispersion coefficient with frequency for various values of particle 


relaxation time. 


